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We formulate an adiabatic approximation for the imaginary-time Schrodinger equation. The ob¬ 
tained adiabatic condition consists of two inequalities, one of which coincides with the conventional 
adiabatic condition for the real-time Schrodinger equation, but the other does not. We apply this adi¬ 
abatic approximation to the analysis of Markovian dynamics of the classical Ising model, which can 
be formulated as the imaginary-time Schrodinger equation, to obtain an asymptotic formula for the 
probability that the system reaches the ground state in the limit of a long annealing time in simulated 
annealing. Using this form, we amend the theory of Somma, Batista, and Ortiz for a convergence 
condition for simulated annealing. 


1. Introduction 

An optimization problem is a problem of finding an element of some set that minimizes a real¬ 
valued function called the cost function. In this paper, we consider an optimization problem with 
discrete variables, which is known as a combinatorial optimization problem. The cost function of a 
combinatorial optimization problem is identified with the Hamiltonian of the classical Ising model 
whose ground state is the global minimum. Solving combinatorial optimization problems is difficult 
in general because of the exponential increase of the number of elements with the problem size and 
frustrations in the problem. It is generally very difficult to find the exact solution within a practical 
time. We thus devise algorithms that give an approximate solution. Simulated annealing 1 2 and quan¬ 
tum annealing 3 El 5 are among such approximate algorithms. 

The basic idea of these algorithms is to use a physical process to escape local minima of the cost 
function so that the state approaches the global minimum. In simulated annealing, we introduce a 
time-dependent temperature T(t) as the control parameter. We initially set the temperature to a high 
value and reduce T{t) slowly toward zero, and the system finally reaches the zero-temperature equi¬ 
librium, the ground state that corresponds to the solution of the combinatorial optimization problem. 
Quantum annealing was proposed in an analogy with simulated annealing . 3 In quantum annealing, we 
introduce a time-dependent external magnetic field which induces quantum fluctuations. We reduce 
the external magnetic field from a very large value to zero, similar to simulated annealing in which 
we reduce the temperature. A similar idea, adiabatic quantum computation^ is often used in the liter¬ 
ature of quantum information theory. However, there is a small difference between adiabatic quantum 
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computation and quantum annealing in that adiabatic quantum computation only uses adiabatic time 
evolution, but nonadiabatic time evolution is also considered in quantum annealing. 7 In this paper, we 
consider quantum annealing following adiabatic time evolution, i.e., adiabatic quantum computation. 

The classical-to-quantum mapping discussed in Refs. [8] 9, 10 allows us to express the thermo¬ 
dynamical properties of classical systems in terms of those of quantum systems in the same spatial 
dimension. Using this mapping, a slow change in the temperature in simulated annealing corresponds 
to a slow change in the Hamiltonian in quantum annealing. Therefore, we can study simulated an¬ 
nealing and quantum annealing from the same perspective. However, the mapped quantum state does 
not follow the real-time Schrodinger equation, and its dynamics is represented as the imaginary-time 
Schrodinger equation. 11 When we consider time-dependent quantities, we need careful analyses ow¬ 
ing to the difference of the dynamics. Somma et al. applied this classical to quantum mapping to 
simulated annealing, and rederived the convergence condition originally proved in Ref. |T2j under the 
ordinary adiabatic condition for the real-time Schrodinger equation. 10 However, the real-time adia¬ 
batic condition does not directly apply to the imaginary-time Schrodinger equation, and their analysis 
should therefore be carefully reexamined. 

In this work, we derive formulas for the adiabatic approximation for the imaginary-time 
Schrodinger equation. This was derived before in Ref. TTfor the norm-conserved but nonlinear case as 
well as in Ref. '14; Our approximation is applicable to the norm-nonconserved and linear case, which 
is related to classical Markovian dynamics. We apply this approximation to simulated annealing and 
obtain an asymptotic formula for the probability that the system reaches the ground state at zero tem¬ 
perature. Using this formula, we rederive the rate of convergence to the ground state discussed by 
Somma et al. 

This paper is organized as follows. In the next section, we derive the adiabatic approximation for 
the imaginary-time Schrodinger equation. Then, in Sect. [3j we review classical-to-quantum mapping 
and rewrite classical Markovian dynamics as the imaginary-time Schrodinger equation. Applying the 
approximation discussed in Sect. [2] to the mapped quantum system, we analyze the probability of 
reaching the ground state in Sect. [4] The convergence condition for simulated annealing is rederived 
from the imaginary-time adiabatic condition in Sect. [5] The final section is devoted to the conclusion. 

2. Imaginary-Time Schrodinger Equation and Its Adiabatic Approximation 

Let us consider the imaginary-time Schrodinger equation 

- i m) = mm))- w 

at 

We consider the time development of a system following this equation in the time scale r, 0 < t < t. 
We scale the time as s = f/r, where s starts from 0 and ends at s = 1. Then, Eq. dTJ reads 

- = tH(s)\i//(s)). (2) 

as 

Note that the norm of the wave function is not conserved, and (tj/(s)\i//{s)) depends on s. 
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Following Refs. [13] and [141 we expand the wave function in terms of the set of instantaneous 
eigenstates, H(s)\j(s)) = Ej(s)\j(s)), as 

l<AO)> = ^Cj(s)\j(s)) = ^ e~ T ^ U) Cj(s)\j(s)}, (3) 


where the second equality defines cjs) with 


4>j(s) = f ds'Ej(s'). 
Jo 


(4) 


We assume Eo(s) = 0 with an application in later sections in mind. The imaginary-time Schrodinger 
equation © can be rewritten in terms of the coefficients as 


ds 


k*j 


d IA s l _ V 


A jk(s) 




(5) 


where A jk(s) = Ej(s)-EJs). Integration of this differential equation and multiplication of the resulting 
expression by g _T W) yield 


WO , V f AJ r T*,«) <W\ dH dV\ k ^y> r ,s 


Cj(s) = c ;(0)e~ 


+ e 




ds e 


A/*(^) 


c k (s ). 


( 6 ) 


Let us solve this integral equation iteratively, i. e ., an asymptotic expansion for very large r. The 
initial condition is that co(0) for the ground state is of (9(1), and the other coefficients are much smaller 
or even zero. Then, the zeroth-order solution Cj, which is obtained by ignoring the integral part in 
Eq. ©, is 

4%) = c 0 (0), cf^ = Cj ( 0)e~ T ^ s \ (7) 

Insertion of these relations into Eq. © gives 


c% Q) (s) = Cj( 0)e~ T W s) + e~ T ^ s) ^ c*(0) f ds' 

k*j 


'n {j(s')\^Jp\k(s')) 


Ajk(s') 


= c 0 (0)e- T ^) 

Integration by parts leads to 

„(l) 


f ds', 

Jo 


,T 0 j(s ) 




Ajo(s') 


+ 0{e T ). 


( 8 ) 


(9) 


c^ 0) (s) = co(0)e~ T ^ 


Sajis') 


,,O'(^)I^P|0(^)> 


-- Cds'e^A 
r Jo ds' 

c 0 (0) (M\^jr\0(s)) 


A/o(s') 2 


OV)I^IO(s')) 


A/o(s') 2 


T A j0 (s) 2 


+ 0{t z ) 


c 0 (0)A ,(s) 


+ 0(r~ 2 ). 


( 10 ) 

( 11 ) 

( 12 ) 
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CnV) = C 0(0) - 


co(0) 


5f 


els' 




Am(s ') 3 


+ 0{t 2 ) 


(13) 


= c 0 (0) |l - ^ ^ ds' B(s')J + 0(r~ 2 ). (14) 

Equations (fT2T) and (fl4l) represent the adiabatic approximation for the imaginary-time Schrddingcr 
equation. 

Equations (fl2l) and (fl4l) suggest that the adiabatic condition for the imaginary-time Schrddingcr 
equation is 


Aj{s) 


« 1, 


fo ds'B(s') 

T 


« 1, 


(15) 


the former of which coincides with the conventional adiabatic condition of the real-time Schrddingcr 
equation. 15 We must be careful, however, that the norm of the wave function is not conserved, and 
hence |c ; (s)| 2 does not directly represent the probability. We shall come back to this point later. 


3. Master Equation Expressed as the Imaginary-Time Schrodinger Equation 


Nonequilibrium dynamics of the Ising model following the master equation can be rewritten as 
the imaginary-time Schrddingcr equation as described in Refs. I Ojand! IJ The master equation is 

1 dP lT (s) 
t ds 


= J] W IT( As)Pa>{s), 


(16) 


where we have scaled the time as s = t/r as before, cr is a set of N Ising spins cr = {o"i, o" 2 ,..., cr ; v), 
and Pcr(s) is the probability that the system is in state cr at scaled time s. We have the Ising model 
with the Hamiltonian Hq(ct) in mind, which is reflected in the transition matrix W aa -’(s) implicitly. 
Note that the transition matrix W (nr ds) may be time-dependent through the time dependence of the 
temperature T(s) or its inverse f3(s). 

Suppose that the transition matrix follows the detailed balance condition 

-P(s)Ho(.cr) j 


W (7cr ds)P^(s) = W^MP^is) 


P ( °\s) = - 


(17) 


The right eigenvalues of the transition matrix are denoted as To(= 0) > Ai > X 2 > ■ ■ ■. The lead¬ 
ing eigenvalue/eigenvector corresponds to thermal equilibrium, which does not change with time as 
suggested by A 0 = 0. 

The following ‘similarity transformation’ is the key to mapping the classical nonequilibrium dy¬ 
namics to quantum mechanics E3]Di)JD il 


H S a(s) = -e^ s)Ho W{s)e-^ s)H ° 

(18) 

As)) = e^ (s)Ho ^ P (T \cr), 

(19) 


cr 


where W(s) is a 2 V x 2 ,v matrix with elements W ITrT '(s), and Hs,\(s) is also a matrix. Note that |t/4.v)) 
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is not normalized. It is easy to see that this H$ A (s) is Hermitian, and can therefore be regarded as 
a quantum-mechanical Hamiltonian. Two matrices, W(s) and H$ A (s), share the spectrum and eigen¬ 
states, up to a trivial factor or sign, 

W(s)|T„(T» = A n (s)\A n (s)) (20) 

H SA (s)\cf> in \s)) = = -A n (s)\cf> {n \s)) (21) 

|0 (m) (a)> = e^ H °\A n (s)), (22) 

as can be verified from Eqs. (flSl ) and (fl9l ). The vectors |i//(.v)> and \(p (n \s)) are not normalized. The 
normalized eigenvector of H$ A (s) will be denoted as |«sa(s)>- I n particular, the normalized ground 


state is 


e ~\fKs)H 0 _ 

|0sa(t>> = — ~j=— |o->, 


(23) 


which corresponds to thermal equilibrium having Aq(s) = 0 and consequently E^(.v) = 0. The expec¬ 
tation value of an arbitrary matrix diagonal in the cr-basis by the ground state of H$ A (s) is equal to the 
expectation value by the Boltzmann distribution. This suggests that thermal fluctuations are mapped 
to quantum fluctuations of the ground state. 

From the master equation (fT6l) . |i/d.v)> can be verified to satisfy the following differential equation, 


- ^#(s)> = t|// S a(.v) - j |tA(T)>, 


(24) 


where ji is for d/3/ds. This is a type of imaginary-time Schrodinger equation with the effective Hamil¬ 
tonian 


HUs) = H sa (s) - 

2 T 


(25) 


The normalized instantaneous eigenstate of the effective Hamiltonian will be written as 

H to t(s)\jtot(s)) = E tot (s)\j {ot (s)). (26) 

4. Probability of Reaching the Ground State 


Let us write the spin configuration of the ground state of Ho as |cr G ). The probability that the 
system reaches the ground state at time s is 


Po-q (s) = <CT G | ^ P o-(A')|cr> 

cr 

= (cTG\e-^ s)Ho \Hs)) 
- e~'^ (s)EG {cT G \ifj(s)). 


(27) 

(28) 
(29) 


This expression can be decomposed as 


P<t q (s) = e ^(cr G \ks A {s))(k SA (s)\j {ot (s))(j tot (s)\il/(s)). 


(30) 
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Now, we assume that the temperature is controlled such that it reaches T — 0 ((3 —* oo) at s = I and 
that the ground-state energy of Ho is also zero, Eq = 0. Then, the instantaneous eigenstate of H$a(s) 
at s = 1 is the ground state, so (ct g |& S a( 1)) = <N,o- We therefore have 


V G (1) = 2<°sA( 1 )l;tot(i)X;tot(i)l<A(i)>. 


(31) 


According to the definition (1251) and perturbation theory, the instantaneous eigenstate of the total 
Hamiltonian is related in the large-r limit to H$a as 


ly'tot} - Usa) - ^ ^ |/sa) 




(IsaIHoUsa) 

~ E^(s) 


+ 0{t 2 ). 


"SA V 


We thus have 


and 


<0 SA (j)ljtot(s)> = - 


<0 S A(s)|0tot(j)> = 1 +0(T 2 ), 
$ <Osa|//oIjsa) 


+ 0(t~ 2 ) u* 0). 


E^lis) - Efl(s) 


Then, from Eq. (I3T1) . 


7V G (i) = <0 tot (i)kKi)> - f 2 


P v 1 <Osa|//oI./sa> 


■(jumw)+ o(t~ 2 ). 


(32) 


(33) 


(34) 


(35) 


2r U - 4>> 

The asymptotic expansions of Eqs. (fT21) and (fl4l) developed in Sect. [2] for the imaginary-time 
Schrodinger equation tell us that 


<0 tot (.s-)|<A(.s-)> = c‘ ot (0)|l -1 




B tot (s')ds' | + 0(t 2 ) 

T Jo 

0 w “ tot + 0 ( t ~ 2 ), 


from which we have 


where 


/V G (l) = c™(0)|l 


If 1 

T JO 


B tot (s)ds I +0(t 2 ), 


BUs) — 


(jtot(^)l% M |0tot(^)> 


Vo (E^(s)~E[ 0 ° t V)) 3 

Since the difference between H tot and Hsa is of (9(r _1 ), we finally obtain 


P aa ( l) = cf(0)|l 


if 1 

T Jo 


B S A(s)ds] + 0 (t 2 ), 


where 


Bsa(s) = 2 
j* 0 


Osa(s)I%^|0 S a(s)> 


£ s>) 3 


(36) 

(37) 

(38) 

(39) 

(40) 

(41) 












J. Phys. Soc. Jpn. 


Full Paper 


and we have replaced Cg 0t (l) by Cg(l) because the difference of these coefficients is of <9(r -2 ) accord¬ 
ing to Eq. (1331) . We hereafter assume Cq A (0) = 1, which indicates that the initial state was the ground 
state of //sa( 0 ), i.e., the thermal equilibrium state at the inverse temperature /?(0) . 


We have also taken into account the fact that the ground-state energy of H$a is £i°i(s) = 0, 


SA V 


H sa (s) 


e -■ 




Z |cr) 


= 0 . 


(42) 


In order to simplify the expression for Z?saG), we take the derivative of the above equation with respect 


to s. 


dH SA (s ) 


ds 




Z |cr) 

CT J 


\ ( 
= H sa (s) 




The projection of this equation to |./sa(t)) gives 


/ • ( t \\ *"SA‘ 

OsaG)I— j-|0 saG)> = 




ds 

from which we have the simplified expression 

0 2 


(jsA(s)\Ho\0 SA (s)), 


Bsa(s) = -r Z 


I<7sa(-v)|2 1 /o|Osa(-v))I 2 

Eg(s) 


(43) 


(44) 


(45) 


j* o 

5. Convergence Condition of Simulated Annealing 

We are now ready to analyze the problems of the analysis in Somma et al.W- They used the 
classical-to-quantum mapping described in Sect. [3] to rewrite classical nonequilibrium dynamics 
as quantum mechanics. Then they applied the conventional adiabatic condition for the real-time 
Schrodinger equation to derive a differential equation for the temperature variable of the original 
classical system. By solving this differential equation, they ‘rederived’ the Geman-Geman 12 condition 

pN 


T(t) 


log t’ 


(46) 


for the original classical dynamics of the Ising model to reach the ground state with probability close 
to unity in the limit of the long time scale, t » <9(1/A), where A is |/li(T)| = E^(.v) in our notation. 
The quantity in the numerator p is an (9(1) constant. 

There are two points of incompleteness in their argument. First, we have to use the imaginary- 
time Schrodinger equation to analyze classical dynamics, not the real-time Schrodinger equation. The 
adiabatic conditions of these two cases have subtle differences as discussed in detail in Sect. |2] The 
second problem is that they did not use the exact expression for the mapped Hamiltonian H saG) 
defined in Eq. (ITSl) but replaced it by a simpler form with the coefficient of the transverse-field term 
being constant, 


H' q = H 


Z' 


A’ 


(47) 


-x/ 

j 

in their notation, where ^ 

Let us discuss the second point first since it is not a very serious one. According to Eq. (fl8l) . the 
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mapped Hamiltonian H$a is a generalized transverse-field Ising model where the coefficient of the 
transverse field generally depends on the spin configuration. For example, the simplest case of the 
one-dimensional Ising model is mapped t(P 

N 1 TV j TV 

Hsa = j -jtanhW^^ + , ~ 2 cosh 28/ Z (cosh ^ 7 ~ sinh 2 #/^.,,^,)^. (48) 

;'=i p 7=i 

(/""-dependence exists in the coefficient of erh Nevertheless, for the purpose of evaluation of the small¬ 
est energy gap, it is allowed to replace the coefficients by their smallest values, which depend on /3 
exponentially as e~ p @(= x)- The reason for this is that the evaluation of the smallest energy gap using 
the Hopf theorem, 16 as discussed in Somma et all — and as described in detail in Lemma 3.3 of Morita 
and Nishimori, 13 uses only the smallest values of the off-diagonal elements. Thus, the resulting general 
lower bound of the energy gap 

A (s) = 4a<t) > a ^/Ne~ 2(pf}(s)+C)N , (49) 


where a and c are ^-independent positive constants, can be used in the present context. 

The first point regarding the difference between imaginary-time and real-time Schrodinger dy¬ 
namics must be taken more seriously, for which reason we have developed a theory of the previous 
sections. If we are allowed to ignore higher-order terms than the first order in r _1 , which itself needs 
verification rigorously speaking, the condition that the ground-state probability is sufficiently close to 
unity is, according to Eqs. (1401) and (1451) . 


If 

T Jo 


0 2 


Bsa(s) = — 2 


B$A(s)ds <s 1 

IOsa(s)|//q|0sa(s))| 2 
77 ( 7 ) / 


(50) 

(51) 


j *0 

To satisfy this condition, the largest term in the above sum (with j = 1) must be very small. If we 
replace the denominator of the expression for Z?saG) by its smallest value in Eq. (l49l) and the matrix 
element in the numerator by its upper bound, a constant times the system size pN, we obtain the 
following condition: 


fo5) 2 ^* = «« 1, (52) 

aViv Jo 

where we have restored the original time scale t = st. The dot over f3 now denotes the derivative with 
respect to t. We next take the limit of the infinite time scale, r —» 00 , which is the situation for which 
the Geman-Geman condition was originally derived. Then, only the upper bound of the above integral 
relation is changed to infinity provided that p is a function of t only, /. e ., without r-dependence. For 
the resulting condition 


4 e 2cN p 2 N 2 


r*oo 

I 0) 2 e 2ppN dt = 6 « 1 

Jo 


(53) 


flViv 

to hold, the integrand should approach zero sufficiently quickly in the large-7 limit. More explicitly, 








J. Phys. Soc. Jpn. 


Full Paper 


/3(t) is expected to asymptotically satisfy the differential equation 

4 e 2cN pW 2e2f}pN = b2rl _ e > 

aVN 

with sufficiently small b. By rewriting the above as 

2e cN pN dp 




dt 


e PpN = &r (l«)/2 j 


we solve it for fi(t) as 


2eCN .^PN = J*_p-e)/2 +( S t 

1 - e 


or 




1 / 2 /? 
ppN = -cN + - log(a VfV) - log 2 + log ( -- t il ~ e)/2 + c 


1 - e 


If we keep only the leading-order term for large t, 




(54) 


(55) 


(56) 


(57) 


(58) 


This agrees with Somma et al. except for a small correction e(> 0). Notice that then - pN is our 2 pN. 


6. Conclusion 

We have established adiabatic-theorem-like relations for the imaginary-time Schrodinger dynam¬ 
ics. This was done before in Ref. IT3lfor norm-conserved dynamics, which is not necessarily suitable 
for the analysis of the master equation of classical Markovian dynamics. De Grandi et al . 14 also dis¬ 
cussed this problem. We developed their calculations further to obtain a more compact expression, as 
seen in Eqs. (fT2l) and (fl4l) . The result was applied to studying the validity of the analysis in Ref. QXS 
which rederived the convergence condition of simulated annealing to the target ground state. We have 
found that the conclusion of Ref.[l0[is correct, but the process to reach it needs more careful analyses 
as developed here. Our theoretical framework may also be used to shed new light on the analysis of 
finite-temperature slow dynamics of classical Ising models, e.g., spin glasses. 
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